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GENERATING A SIMULATED FLUID FLOW 
OVER AN AIRCRAFT SURFACE USING 
ANISOTROPIC DIFFUSION 

STATEMENT REGARDING FEDERALLY 5 

SPONSORED RESEARCH OR DEVELOPMENT 

This invention was made with Government support under 
contract NNL08AA08C awarded by NASA. The Govern- 
ment has certain rights in the invention. to 

BACKGROUND 

1. Field 

This application relates generally to simulating a fluid flow 15 
over a computer-generated aircraft surface and, more specifi- 
cally, to estimating a viscous fluid property distribution over 
a computer-generated aircraft surface using an anisotropic 
diffusion technique. 

2. Description of the Related Art 20 

Aerodynamic analysis of an aircraft moving through a fluid 

typically requires an accurate prediction of the properties of 
the fluid surrounding the aircraft. Accurate aerodynamic 
analysis is particularly important when designing aircraft 
surfaces, such as the surface of a wing or control surface. 25 
Typically, the outer surface of a portion of the aircraft, such as 
the surface of a wing, is modeled, either physically or by 
computer model, so that a simulation of the fluid flow can be 
performed and properties of the simulated fluid flow can be 
measured. Fluid-flow properties are used to predict the char- 30 
acteristics of the wing including lift, drag, boundary-layer 
velocity profiles, and pressure distribution. The flow proper- 
ties may also be used to map laminar and turbulent flow 
regions near the surface of the wing and to predict the forma- 
tion of shock waves in transonic and supersonic flow. 35 

A computer-generated simulation can be performed on a 
computer-generated aircraft surface to simulate the fluid 
dynamics of a surrounding fluid flow. The geometry of the 
computer-generated aircraft surface is relatively easy to 
change and allows for optimization through design iteration 40 
or analysis of multiple design alternatives. A computer-gen- 
erated simulation can also be used to study situations that may 
be difficult to reproduce using a physical model, such as 
supersonic flight conditions. A computer-generated simula- 
tion also allows a designer to measure or predict fluid-flow 45 
properties at virtually any point in the model by direct query, 
without the difficulties associated with physical instrumenta- 
tion or data acquisition techniques. 

In some cases, a computer-generated simulation includes a 
computational fluid dynamics (CFD) simulation module used 50 
to predict the properties of the fluid flow. A CFD simulation 
module estimates the properties of a simulated fluid flow by 
applying an algorithm or field equation that estimates the 
interaction between small simulated fluid volumes, also 
referred to as fluid cells. Because a single CFD simulation 55 
module may include millions of individual fluid cells, the 
complexity of the relationship between fluid cells can have a 
large effect on the computational efficiency of the simulation. 
Complex CFD simulation modules can be computationally 
expensive and require hours or even days to execute using 60 
high-perfonnance computer processing hardware. 

To reduce the computational burden, in some instances it is 
desirable to use a CFD simulation module that simplifies the 
fluid dynamics and produces a fluid simulation that can be 
solved more rapidly. For example, for fluid flows that are 65 
relatively uniform or are located away from an aircraft sur- 
face, a simplified simulation that minimizes or ignores certain 
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fluid dynamic phenomena can be used. In the examples dis- 
cussed below, a simplified simulation may ignore fluid 
dynamic contributions due to fluid viscosity, which, in some 
cases, have little effect on the overall behavior of the fluid. By 
using a simplified simulation, processing time may be 
improved. 

In other situations, where the fluid flow is not as uniform, it 
may be necessary to use a CFD simulation module that is 
more sophisticated and capable of accurately predicting the 
fluid properties, using more complex fluid dynamics. In the 
examples discussed below, a more sophisticated simulation 
may account for dynamic contributions due to fluid viscosity. 
Under certain conditions, the more sophisticated simulation 
may be more accurate, particularly for portions of the fluid 
flow where fluid viscosity affects the results. However, more 
sophisticated simulation modules are also likely to require 
more computing resources, and therefore may require more 
time to solve. 

SUMMARY 

One exemplary embodiment includes a computer-imple- 
mented method of generating a fluid-flow simulation over a 
computer-generated aircraft surface. The surface is com- 
prised of a surface mesh of polygons. A boundary-layer fluid 
property is obtained for a subset of the polygons of the surface 
mesh of polygons. The subset represents a portion of the 
computer-generated aircraft surface. A pressure-gradient 
vector is determined for a selected polygon, the selected 
polygon belonging to the surface mesh of polygons but not 
one of the subset of polygons for which a boundary -layer fluid 
property has been obtained. A minimum diffusion rate is 
defined along the direction of the pressure gradient vector of 
the selected polygon. A maximum diffusion rate is also 
defined in a direction perpendicular to the direction of the 
pressure gradient vector of the selected polygon. A diffusion- 
path vector is defined between a point in the selected polygon 
and a neighboring point in a neighboring polygon, the neigh- 
boring polygon belonging to the subset of polygons for which 
a boundary-layer fluid property has been obtained. An 
updated fluid property is determined for the selected polygon 
using a variable diffusion rate, the variable diffusion rate 
based on the minimum diffusion rate, maximum diffusion 
rate, and angular difference between the diffusion-path vector 
and the pressure gradient vector. 

DESCRIPTION OF THE FIGURES 

FIG. 1 depicts a computer-generated fluid flow applied to a 
computer-generated aircraft surface. 

FIG. 2 depicts an exemplary process for propagating a 
boundary-layer fluid property over a computer-generated air- 
craft surface using anisotropic diffusion. 

FIG. 3 depicts an exemplary surface mesh and fluid-flow 
mesh for simulating a fluid flow over a computer-generated 
wing surface. 

FIG. 4 depicts a surface mesh of polygons representing a 
wing surface and strip areas representing a portion of wing 
surface. 

FIG. 5 depicts results of a boundary-layer CFD simulation 
performed only for strip areas of a wing surface. 

FIG. 6 depicts exemplary diffusion paths on a surface mesh 
of polygons representing a computer-generated wing surface. 

FIGS, la and lb depict an exemplary fluid flow around a 
wing surface. 
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FIG. 8 depicts an exemplary pressure distribution over a 
computer-generated wing surface using isotropic diffusion 
and a finite-volume solver. 

FIG. 9 depicts an exemplary pressure distribution over a 
computer-generated wing surface using an inviscid CFD 5 
simulation module. 

FIG. 10 depicts exemplary pressure distribution over a 
computer-generated wing surface using anisotropic diffusion 
and a finite-volume solver. 

FIG. 11 depicts an exemplary computer system for simu- to 
lating a fluid flow over a computer-generated aircraft surface. 

FIG. 12 depicts an exemplary computer network. 

The figures depict one embodiment of the present inven- 
tion for purposes of illustration only. One skilled in the art will 
readily recognize from the following discussion that alterna- 15 
tive embodiments of the structures and methods illustrated 
herein can be employed without departing from the principles 
of the invention described herein. 

DETAILED DESCRIPTION 20 

As discussed above, a computer-generated simulation can 
be used to analyze the aerodynamic performance of a pro- 
posed aircraft surface, such as a wing or control surface. 
Using known geometry modeling techniques, a computer- 25 
generated aircraft surface that represents the outside surface 
of the proposed aircraft can be constructed. FIG. 1 depicts an 
exemplary computer-generated aircraft surface of the Space 
Shuttle orbiter vehicle, external tank, and twin solid-fuel 
rocket boosters. A CFD simulation module has been applied 30 
using the computer-generated aircraft surface of the Space 
Shuttle orbiter to predict the fluid properties of an exemplary 
fluid flow. 

As shown in FIG. 1, the results of the simulation can be 
visually represented as shaded regions on the computer-gen- 35 
erated aircraft surface of the Space Shuttle. Different shades 
represent the predicted pressure distribution resulting from 
the simulated fluid flow. In FIG. 1, transitions between the 
shaded regions represent locations of predicted pressure 
change across the surface of the Space Shuttle. Similarly, 40 
different pressures in the surrounding fluid flow are repre- 
sented as differently shaded regions. 

In FIG. 1, the simulation of the fluid flow is visualized by 
depicting the predicted pressure distribution. However, the 
simulation may be visualized using other fluid properties, 45 
including surface velocity, air temperature, air density, and 
others. Additionally, the simulation may be used to visualize 
locations of developing shock waves or transitions between 
laminar and turbulent flow. 

The simulation allows the designer or engineer to evaluate 50 
the performance of the aircraft geometry for various flow 
conditions. If necessary, changes can be made to the aircraft 
geometry to optimize performance or eliminate an unwanted 
aerodynamic characteristic. Another simulation can be per- 
formed using the modified geometry, and the results can be 55 
compared. To allow for multiple design iterations, it is advan- 
tageous to perform multiple simulations in a short amount of 
time. However, as described above, there is a trade-off 
between speed and accuracy of the simulation depending on 
the type of CFD simulation module used. 60 

As mentioned above, a computer-generated simulation 
typically represents a fluid flow as a three-dimensional fluid- 
flow mesh of small volumes of fluid called fluid cells. The 
shape of the fluid cells can vary depending on the type of 
fluid-flow mesh used. A CFD simulation module predicts the 65 
fluidic interactions between the fluid cells in the fluid-flow 
mesh, using a fundamental algorithm or field equation. 
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The speed and accuracy of a CFD simulation module 
depends, in part, on the field equation used to predict the 
interaction between the fluid cells. In some instances, the field 
equation simplifies the relationship between fluid cells by 
ignoring or minimizing certain dynamic contributions. These 
field equations are typically less complex, and therefore are 
more computationally efficient. For instance, a simplified 
algorithm called the Euler method may be used to simulate a 
fluid flow when viscous effects can be minimized or ignored. 
Viscous effects of a fluid can be ignored when, for example, 
there is not a significant velocity difference between adjacent 
fluid cells, and therefore shear forces due to internal friction 
or viscosity are minimal. A CFD simulation module that 
ignores or minimizes effects of fluid viscosity can also be 
referred to as an inviscid simulation. 

In other instances, a more complex field equation is used to 
more accurately predict the interaction between the fluid 
cells. For example, a Navier-Stokes method can be used to 
simulate the pressure and shear forces on the fluid cells. 
Unlike the Euler method mentioned above, the Navier-Stokes 
method accounts for the effects of viscosity and offers a more 
accurate simulation of a fluid flow. A simulated fluid flow that 
accounts for effects due to fluid viscosity can also be referred 
to as a viscous simulation. 

However, the improved accuracy of the Navier-Stokes 
method comes at the cost of increased computational load, 
and therefore the Navier-Stokes method is generally slower to 
compute than an Euler-based algorithm. In addition, a 
Navier-Stokes method may require more (smaller) fluid cells 
to generate acceptable results, further increasing the compu- 
tational load. Thus, selecting the field equation for a CFD 
module often involves a trade-off between speed and accu- 
racy. In practice, designers may use faster Euler-based CFD 
models to evaluate multiple design iterations and then vali- 
date the final design iteration with a more accurate Navier- 
Stokes-based CFD model. However, if the Navier-Stokes 
CFD simulation reveals a design problem, the entire process 
must be repeated, wasting valuable time and computing 
resources. 

The techniques described below are computer-generated 
simulations that use both viscous and inviscid field equations 
to achieve acceptable accuracy without requiring the compu- 
tational burden of a full viscous CFD simulation (e.g., a 
Navier-Stokes-based CFD simulation). In many simulations, 
there is a region of the fluid flow that can be accurately 
predicted without taking viscous contributions into account. 
For example, regions of the fluid flow that are located away 
from an aircraft or wing surface have a relatively uniform 
velocity profile. Therefore, an inviscid simulation using, for 
example, an Euler-based analysis, can be used to accurately 
predict the behavior of the fluid-flow region. In other fluid- 
flow regions where there is a less uniform velocity profile, a 
more complex, viscous simulation can be used. 

In fluid-flow regions where a viscous simulation is desired, 
further computational gains can be realized by performing the 
viscous simulation for only a portion of the fluid-flow region. 
For example, a portion of a boundary-layer fluid-flow region 
(near a wing surface) may be represented by one or more strip 
areas. See, for example FIG. 4 , depicting strip areas 404 
extending along the chord direction of the wing and spaced at 
intervals along the span of the wing. In FIG. 4 , the strip areas 
404 are spaced at regular intervals; however, it is not neces- 
sary that they be so. The viscous simulation is performed only 
for the strip-area portions of the boundary-layer fluid-flow 
region, thereby reducing the number of simulation calcula- 
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tions. As a result, the boundary-layer fluid properties are only 
known for strip-area portions of the boundary-layer fluid- 
flow region. 

Using the techniques described below, the predicted 
boundary-layer fluid properties for a portion of a boundary- 5 
layer region (represented by, for example, strip areas 404) can 
be extended to the other portions of the boundary-layer 
region. Specifically, a diffusion technique can be used to 
propagate a predicted boundary-layer fluid property for the 
strip areas 404 to other areas of the computer-generated air- 10 
craft surface, between the strip areas 404. In this way, the 
techniques described herein may provide an accurate predic- 
tion of the fluid properties of an entire fluid-flow region 
without having to perform a CFD simulation for the entire 15 
boundary-layer fluid-flow region. 

1 . Exemplary Process for Propagating a Fluid Property Using 
Diffusion Techniques 

FIG. 2 depicts an exemplary process 200 for propagating 
the results of a boundary-layer CFD simulation module for a 20 
computer-generated aircraft surface. The results of the 
boundary-layer CFD simulation module include one or more 
boundary-layer fluid properties (e.g., pressure, temperature, 
surface velocity, boundary-layer thickness). For purposes of 
this discussion, a computer-generated wing surface is used. 25 
However, the technique can be applied to nearly any type of 
computer-generated aircraft surface, including a computer- 
generated surface for an aircraft fuselage, nacelle, shroud, or 
other structure. 

The exemplary process 200 uses a diffusion technique to 30 
propagate a fluid property over the computer-generated air- 
craft surface. The term diffusion, as it is used below, relates to 
a technique for estimating fluid property values over a com- 
puter-generated surface using fluid property values that are 
(initially) known for only a portion of the computer-generated 35 
surface. The estimation is based, in part, on a modified ver- 
sion of the heat equation, which is also referred to herein as a 
diffusion equation. Traditionally, a diffusion equation can be 
used to predict heat diffusion or conduction through a mass or 
body of material. However, rather than predict the physical 40 
diffusion of heat energy, the exemplary process 200 uses a 
diffusion equation to propagate or distribute fluid property 
values from one portion of a computer-generated surface 
(where the fluid properties are known) to other portions of the 
computer-generated surface (where the fluid properties are 4 - 
not known). 

An exemplary two-dimensional, isotropic diffusion equa- 
tion that can be used to diffuse a fluid property u over a 
distance x and y can be expressed as: 


6u _(S 2 u 6 2 u 1 Equation 1 

5t 4 i).v- + 6y 2 ]' 

55 

where a is a positive constant representing the diffusion rate 
and t is time. The diffusion rate a (also called a conduction 
constant) determines the rate of diffusion over a given dis- 
tance, also called a diffusion path. As applied in the exem- 
plary process 200 below, a diffusion path is typically defined 60 
between two points on a computer-generated aircraft surface: 
a first point having an unknown fluid property value and a 
second point with a known fluid property value. In general, a 
large diffusion rate tends to produce a more uniform fluid- 
property distribution between the two points along a diffusion 65 
path. That is, the unknown fluid property value at the first 
point will be closer to the known fluid property value at the 
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second point. Conversely, a small diffusion rate tends to pro- 
duce a less uniform fluid-property distribution between two 
points along a diffusion path. 

Typically, process 200 is performed as part of a computer- 
generated simulation for a computer-generated aircraft sur- 
face. As mentioned above and discussed in more detail below 
with respect to FIGS, la and lb, a computer-generated simu- 
lation can be performed using two CFD simulation modules. 
In the exemplary process 200, an inviscid CFD simulation 
module is used to predict the fluid properties for an inviscid 
fluid-flow region and a viscous or boundary-layer CFD simu- 
lation module is used to predict the fluid properties of a 
viscous fluid-flow region. 

As shown in FIG. 3, the inviscid fluid-flow region may be 
discretized using a fluid-flow mesh 302. In the following 
example, fluid-flow mesh 302 represents the inviscid fluid- 
flow region using a three-dimensional array of fluid cells, 
each fluid cell associated with one or more inviscid-fluid 
properties. As explained in more detail below, an inviscid 
CFD simulation module (using, for example, field equation 
13, below) can be used to determine inviscid-fluid properties 
for each fluid cell of the fluid-flow mesh 302. 

FIG. 3 also depicts a computer-generated aircraft surface, 
specifically a wing surface, discretized using a surface mesh 
of polygons 304. While FIG. 3 depicts a surface mesh of 
triangular polygons, other types of polygons can also be used. 
FIG. 3 also depicts an exemplary strip area 404 and subset of 
polygons 420, discussed in more detail below with respect to 
FIGS. 4, 5, and 6. 

The boundary -layer fluid-flow region is discretized using a 
two-dimensional array of points within the surface mesh of 
polygons 304. As shown in FIG. 3, the boundary -layer region 
may be represented using boundary-layer prediction points 
306 located on the computer-generated surface. In this 
example, each boundary-layer prediction point 306 corre- 
sponds to a polygon in the surface mesh of polygons 304. The 
boundary-layer prediction points 306 can be located virtually 
anywhere with respect to a corresponding polygon. However, 
for ease of explanation, the boundary-layer prediction points 
306 are depicted as being located near the centroid of each 
corresponding polygon. While the computer-generated wing 
surface is a three-dimensional representation, the surface 
mesh of polygons 304 and corresponding boundary-layer 
prediction points 306 can be treated as a two-dimensional 
space projected onto the three-dimensional computer-gener- 
ated wing surface. That is, each polygon or boundary-layer 
prediction point 306 can be identified and located with 
respect to the computer-generated wing surface using two 
coordinate values (x and y). 

Discussed in more detail below with respect to equation 2, 
FIG. 3 also depicts exemplary selected polygon 630 and a set 
of fluid cells 320, close in proximity to selected polygon 630. 

The following operations can be performed as an exem- 
plary process 200 for diffusing a boundary-layer property 
over a computer-generated surface. In operation 202, a 
boundary-layer fluid property for a subset of polygons is 
obtained. As shown in FIG. 4, the subset of polygons 420 may 
be defined as the surface-mesh polygons within a strip area 
404. For simulations using a computer-generated wing sur- 
face, as shown in FIG. 5, the strip areas 404 may extend along 
the chord direction or stream-wise across the computer-gen- 
erated wing surface. Also shown in FIG. 5, multiple strip 
areas 404 may be spaced at a regular distance along the span 
of the wing. In alternative arrangements, the strip areas 404 
can be oriented along different directions and placed at vari- 
ous locations along the computer-generated wing surface. 
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In this example, a strip area 404 is used to define the subset 
of polygons 420. However, other shapes or techniques can be 
used to select the subset of polygons 420, as long as the 
number of polygons in the subset of polygons 420 is fewer 
than the number of polygons in the surface mesh of polygons 5 
304. For example, the subset of polygons 420 may be defined 
as those polygons located within a radial distance of a point 
on the computer-generated aircraft surface. In another 
example, the subset of polygons 420 may be one polygon. 

To reduce the number of simulation calculations, the 10 
boundary-layer simulation is only performed for boundary- 
layer prediction points 306 associated with the subset of 
polygons 420. Using the example shown in FIG. 4, the bound- 
ary-layer simulation is performed for each boundary-layer 15 
prediction point 306 associated with a polygon of the subset 
of polygons 420 (e.g., polygons within strip areas 404). Note, 
a boundary-layer simulation is not performed for boundary- 
layer prediction points 306 associated with polygons that are 
not part of the subset of polygons 420 (e.g., polygons between 20 
strip areas 404). 

As explained in more detail below with respect to FIGS, la 
and lb, a boundary-layer CFD simulation (using, for 
example, field equations 14 and 15 , below) can be used to 
determine a boundary-layer fluid property for each boundary- 25 
layer prediction point 406 associated with a polygon of the 
subset of polygons. 

FIG. 5 depicts exemplary results of operation 202. That is, 
FIG. 5 depicts the results of a boundary- layer CFD simulation 
performed only for a subset of polygons 420 (within strip 3< 1 
areas 404). The shaded areas in FIG. 5 represent the portions 
of the boundary-layer fluid-flow region where the boundary- 
layer fluid properties are known. Changes in the degree of 
shading represent variations in the boundary-layer fluid prop- , _ 
erty values. Note that the grey regions (between strip areas 
404) represent portions of the boundary-layer fluid-flow 
region where the boundary-layer fluid properties are not 
known. 

In operation 204, a pressure gradient vector Vp is deter- 40 
mined for a selected polygon. In the examples shown in FIGS. 

3 and 6, the selected polygon 630 is a polygon of the surface 
mesh of polygons 304 but is not one of the subset of polygons 

420. The direction component of pressure gradient vector Vp 45 
indicates the direction of greatest pressure change with 
respect to the fluid flow near the selected polygon 630. 

The pressure gradient Vp can be determined using the 
inviscid fluid-properties of the fluid-flow mesh 302. As 
explained above, an inviscid CFD simulation module can be 
used to predict the inviscid fluid properties of fluid cells in the 
fluid-flow mesh 302. Typically, the inviscid fluid properties 
are known for at least the fluid cells (of the fluid-flow mesh 
302) that are close in proximity to the computer-generated 
aircraft surface used for the simulation (e.g., the computer- 55 
generated wing surface). In contrast, as mentioned above, the 
boundary-layer fluid properties are (initially) only known for 
the subset of polygons 420 within strip area 404. 

Changes in the inviscid fluid pressure can be indicative of 
changes in the boundary-layer fluid properties. Therefore, an 60 
inviscid pressure gradient vector can be used as a suitable 
proxy for the gradient of the boundary-layer property being 
diffused. 

In one example, the pressure gradient vector Vp is deter- 65 
mined by identifying a set of fluid cells 320 of the fluid-flow 
mesh 302 that are close in proximity to the selected polygon 


630. (See FIG. 3.) Generally, a two-dimensional pressure 
gradient vector in an x-y space can be defined as: 

— . /Sp 6p\ Equation 2 


where p is the pressure of the inviscid fluid flow near the 
surface of the wing (represented by the identified set of fluid 
cells 320 close in proximity to the selected polygon 630). 
Because the wing is a three-dimensional shape, the proximate 
inviscid fluid flow represented by the identified set of fluid 
cells 320 will also be three-dimensional. However, for the 
sake of simplicity, an inviscid fluid pressure distribution cor- 
responding to the identified set of fluid cells 320 can be 
treated as a two-dimensional distribution mapped to the com- 
puter-generated aircraft surface. Because, in many cases, the 
inviscid fluid pressure distribution is different for different 
areas of the computer-generated aircraft surface, a pressure 
gradient should be determined for each selected polygon 
within the surface mesh of polygons. 

As explained in more detail below with respect to opera- 
tions 206 and 212, the pressure gradient vector Vp is used to 
determine the primary and secondary diffusion directions. 

For example, the direction of the pressure gradient vector 

Vp may be used as the first primary direction dj . The direction 

perpendicular to the pressure gradient vector Vp may be used 
as the second primary direction d 2 . 

In operation 206, maximum and minimum diffusion rates 
are defined ( a max , a mi „, respectively). As discussed above, 
the diffusion rate determines, in part, the amount of diffusion 
along a diffusion path. In this example, a minimum diffusion 
rate a„„„ is selected for diffusion paths along the first primary 

direction dj (the direction of the pressure gradient vector Vp). 
A maximum diffusion rate a max is selected for diffusion paths 
along the second primary direction d 2 (the direction perpen- 
dicular to the direction of the pressure gradient vector Vp). In 
some cases, the maximum diffusion rate a max is 1,000 to 
1 ,000,000 times the minimum diffusion rate a mirl . 

In some cases, when the pressure gradient vector Vp is 
small or nearly zero, a single diffusion rate can be used. If the 
diffusion rate is equal or nearly equal in all directions, the 
maximum and minimum diffusion rates (ct max , a min ) may be 
treated as equal, or one of the diffusion rates (a max ) can be 
used as the sole diffusion rate. For example, as discussed in 
more detail below with respect to equation 9, when the pres- 
sure gradient vector Vp is small or nearly zero, the diffusion 
rate does not depend on direction, and isotropic diffusion can 
be applied using the maximum diffusion rate ( o . max ). 

In operation 208, a diffusion path vector cf is defined 
between the selected polygon and a neighboring polygon. 
The neighboring polygon may or may not be immediately 
adjacent to the selected polygon. FIG. 6 depicts an exemplary 
selected polygon 630 and neighboring polygon 620. The dif- 
fusion path vector 640 cf between the selected polygon 630 
and neighboring polygon 620 can be defined as a vector 
between the respective boundary-layer prediction points 306 
of the selected polygon 630 and neighboring polygon 620. 

The direction d di g- use of the diffusion path vector 640 cf can be 
expressed as an angle from horizontal. 
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In FIG. 6, the selected polygon 630 is a polygon outside of 
the strip area 404. Therefore, selected polygon 630 is not one 
of the subset of polygons for which a boundary-layer fluid 
property is known. In this example, neighboring polygon 620 
is inside of the strip area 404, and therefore is one of the subset 5 
of polygons 420 for which a boundary-layer fluid property is 

known. Thus, in this example, the diffusion path vector 640 cf 
defines a path from a portion of the boundary-layer fluid-flow 
region where the fluid properties are known to a portion of the 
boundary layer fluid-flow region where the fluid properties 10 
are not known. 

Other diffusion paths are also possible. For example, FIG. 

6 depicts an alternative selected polygon 631 and alternative 
neighboring polygon 621. Here, both the selected polygon 
631 and neighboring polygon 621 are outside of the strip area 15 

404. The diffusion path vector 641 a defines a path between 
two portions of the boundary layer fluid-flow region where 
the fluid properties are unknown. 

In operation 210, a variable diffusion rate a v is determined. ., (| 
The variable diffusion rate a v depends on the direction d di ^ use 

of the diffusion path vector a , the maximum diffusion rate 
a max , minimum diffusion rate a mi „, and first primary direc- 
tion d : . As explained above, the minimum diffusion rate a min 
applies to diffusion along the first primary direction d, (the 25 

direction of the pressure gradient vector Vp). The maximum 
diffusion rate a max applies to diffusion along a second pri- 
mary direction d 2 (the direction perpendicular to the direction 

of the pressure gradient vector Vp). For diffusion in directions 30 
other than the first or second primary directions (d 15 d 2 ), an 
elliptical relationship can be used to determine the diffusion 
rate associated with a specified diffusion direction d dj)ruse . An 
elliptical relationship for the variable diffusion rate a v can be 
expressed as: 35 

a v =/a mi „ 2 cos 2 (e)+a mm 2 sin 2 (0). Equation 3 

Therefore, the diffusion rate a v is equal to the square root of 
the sum of the square of the minimum diffusion rate a max 
times cosine theta 0 and the square of the maximum diffusion 40 
rate a mUl times sine theta 0 (where theta 0 is the angle between 
the specified diffusion direction d di g use and the first primary 
direction d : ). While equation 3 depicts one exemplary tech- 
nique for determining the variable diffusion rate a . v , other 
techniques can be used to determine the variable diffusion 45 
rate a, using the diffusion path vector d di ^ use , the maximum 
diffusion rate a max , minimum diffusion rate a mirr and first 
primary direction d x . 

In operation 212, a boundary-layer fluid property is deter- 
mined for the selected polygon. A diffusion equation can be 50 
used to determine the boundary-layer fluid property u for the 
selected polygon. An anisotropic diffusion equation can be 
used to account for the changes in diffusion rate due to 
changes in the diffusion direction. An anisotropic diffusion 
equation for fluid property u over a time increment t can be 55 
expressed as: 


du d / d u \ Equation 4 

dt dd\ v dd J’ 60 

where, as shown above in equation 3, a v is a function of the 
maximum diffusion rate a max , the minimum diffusion rate 
a mi „, the first primary direction d x , and the diffusion direction 
d diffuse- Because the first primary direction dj varies depend- 
ing on the direction of the gradient vector Vp, the value of a v 


may also vary for different locations on the computer-gener- 
ated wing surface. Hie diffusion angle d dl ~ t , and the diffu- 
sion distance d may be the direction and magnitude compo- 
nents, respectively, of a diffusion path vector cf. 

As explained in more detail below with respect to dis- 
cretized modeling techniques (Section 2, below), anisotropic 
diffusion may be implemented over a computer-generated 
aircraft surface using two-dimensional coordinates x and y. 
Anisotropic diffusion for a fluid property u over a time incre- 
ment t can also be expressed in terms of two-dimensional 
coordinates x and y: 


du Equation 5 

It = 


d ( 

du ' 

,.i.i 

( du\ 

i. d \ 

( du\ 

d j 

( du\ 

dxV* 

c dx, 

dx ' 



< ayx M 




where a vv , a rv , a VT , and a„, are components of a two-dimen- 
sional diffusion coefficient tensor. The values of a xx , a , a vx , 
and a vv can be determined using the following exemplary 
relationships: 


a xx =a mm COS 2 (<t>): 


Equation 6 


a xy =a yx =a mi „ cos (<}>)sin((j))-a wax cos (<J>)sin(<(>) ; Equation 7 


a yy=CL min cos 2 (<J))+a wa;c sin 2 (4>); Equation 8 

where phi <|> is the angle between the second primary direction 
d 2 and the x-axis. Thus, similar to a v as shown above in 
equation 3, a xx , a vv , a yx , and a yy are also a function of maxi- 
mum diffusion rate a max , the minimum diffusion rate a mjn , 
the first primary direction d 15 and the diffusion direction 

^diffuse' 

In some cases, as mentioned above with respect to equation 

2, a pressure gradient vector Vp that is small or near zero may 
indicate that the diffusion rate does not vary significantly with 
respect to diffusion direction. Thus, when the pressure gradi- 
ent vector Vp is small or near zero, either the maximum or 
minimum diffusion rates ( a max or a mtn , respectively) may be 
used as the diffusion rate a v to predict the fluid property using 
an isotropic diffusion equation. For example, below is an 
exemplary two-dimensional isotropic equation that can be 
used to diffuse a fluid property u over time t along directions 
x and y: 


du d ( du\ d ( du\ Equation 9 

~dt = dx[ av dx) + fry\ Xv ~dy) 


In this special case, the variable diffusion rate a v can be set to 
either the maximum diffusion rate a max , the minimum diffu- 
sion rate a min , or an average of the maximum and minimu m 
diffusion rates. 

Using either equations 4, 5, or 9, a boundary-layer fluid 
property u for the selected polygon can be determined. The 
operations of process 200 can be repeated for multiple 
selected polygons of the surface mesh of polygons. In some 
cases, the operations of process 200 are repeated for each 
polygon (of the surface mesh of polygons 304) that is not one 
of the subset of polygons 420 for which a boundary-layer 
property was obtained using the boundary-layer CFD mod- 
ule. For example, with regard to the examples depicted in 
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FIGS. 3,4,5, and 6, the process 200 can be repeated for every 
polygon in the surface mesh of polygons 304 between strip 
areas 404. 

Equations 4, 5, and 9 are each expressed with respect to a 
small time step dt and are capable of determining a transient 
or time-dependent fluid property value. However, the tran- 
sient or time-dependent results are not necessary to estimate 
the fluid property distribution over the computer-generated 
aircraft surface. In most cases, the steady-state solution is 
used as the estimation of the boundary-layer fluid property u. 
For example, if the differential equations 4, 5, or 9 are solved 
using an explicit method, the explicit method is iterated over 
multiple time steps until a steady-state solution is determined. 
If, instead, the differential equations 4, 5, or 9 are solved using 
an implicit method, a steady-state condition is assumed and 
the boundary-layer fluid property u is solved simultaneously 
for multiple polygons in the surface mesh of polygons. An 
example of explicit and implicit solution techniques is 
described in more detail below in Section 2. 

2. Discretized Modeling Techniques 

The diffusion differential equations 4, 5, and 9 may be 
implemented over a surface mesh of polygons 304 by using a 
numerical approximation technique. As shown in FIGS. 3, 4, 
and 6, each polygon of the surface mesh of polygons 304 may 
have an associated boundary-layer prediction point 306. The 
diffusion equations 4, 5, and 9 can be approximated using an 
explicit finite difference solver that references values of the 
fluid property u at each boundary-layer prediction point 306 
designated by coordinate location i, j. A simplified two-di- 
mensional finite difference solver using central differencing 
for the spatial derivatives can be expressed as: 


= -m"" 1 + Equation 10 

f «i+l, j - 2«, J+M, -I, j K/J-U -2M/J+M/J-1 T“‘ 

H IP + A? J Af + £ ' 

where n is the current time step and n-1 is the previous time 
step, i is the index along the x-direction, j is the index along 
the y-direction, and E is the error. Equation 1 0 can be used for 
a surface mesh of rectangular polygons. A variant of equation 
10 may be used for surface meshes having polygons that are 
not rectangular (e.g., triangular). 

Equation 10 describes an explicit solving method that can 
be used to diffuse a fluid property u among surface-mesh 
points associated with polygons not belonging to the subset of 
polygons 420 (e.g., outside of the strip areas 404). The 
explicit method of equation 10 is typically iterated over mul- 
tiple increments of time t imtil the results converge to a 
steady-state solution. The error term E is included in equation 
10 to account for error due to the discretization of the diffu- 
sion equation. To minimize the error E, the time step incre- 
ment should be sufficiently small. Too large of a time step 
increment will produce inaccurate results and the explicit 
method may become unstable. To maintain stability, the time- 
step increment may be directly proportional to the square of 
the longest dimension of the surface-mesh polygon. 

Because smaller time-step increments require more itera- 
tions before reaching a steady state, the explicit solving 
method using equation 10 may require significant computing 
resources. Although the solution at each time step may be 
relatively easy to calculate, it may take tens of thousands of 
iterations to produce reliable results. 

Alternatively, an implicit method can be used by assuming 
a steady-state condition. For example, at steady-state, the sum 
of the spatial derivatives and the error are zero, allowing: 


« i+1J - 2u iJ + i U-- j+l - 2u iJ + u iJ ~ l n Equation 11 

A^ + A? " °‘ 

5 

Equation 1 1 can be used for a surface mesh of rectangular 
polygons. A variant of equation 1 1 may be used for surface 
meshes having polygons that are not rectangular (e.g., trian- 
gular). Instead of stepping through time, as with the explicit 
to method shown in equation 10, a set of linear equations can be 
constructed. Each linear equation of the set corresponds to a 
boundary-layer prediction point for a polygon of the surface 
mesh of polygons. Mathematically, this can be represented 
as: 

15 

[A ] [u] = 0 , Equation 1 2 

20 where [A] is the matrix of coefficients computed using equa- 
tion 1 1 and [u] is the solution vector that represents the value 
of the diffused fluid property at each boundary-layer predic- 
tion point. For large surface meshes having a large number of 
boundary-layer prediction points, the matrix [A] can become 
25 very large. However, matrix [A] is typically sparse (i.e., popu- 
lated mostly with zeros), and therefore even a large matrix [A] 
may be relatively easy to solve. Using an efficient sparse- 
matrix solver, the implicit method of equation 12 may be 
solved orders of magnitude faster than iterating through the 
30 explicit method shown in equation 10, above. 

While the finite difference techniques exemplified in equa- 
tions 1 0 and 1 2 are relatively straightforward to calculate, the 
results may vary depending on the geometry of the surface- 
mesh polygons. In particular, a surface mesh of polygons that 
35 is less orthogonal than others may not diffuse a fluid property 
equally in all directions. A mesh is said to be more orthogonal 
when there are more lines that can be drawn between the 
centroids of adjoining polygons that intersect the common 
edge between the polygons perpendicularly. 

40 As an alternative, a finite-volume method can be used to 
discretize the diffusion differential equations 4, 5, and 9. In 
one example, a volume cell is defined for polygon of the 
surface mesh. The volume cell may be defined, for example, 
as the geometry of the polygon extruded a nominal thickness. 
45 Finite- volume methods typically provide a more uniform dif- 
fusion because a finite-volume method requires that the 
energy related to the fluid property u be conserved between 
volumes associated with polygons of the surface mesh. 
Therefore, it may be advantageous to use a finite-volume 
50 method instead of finite-difference when discretizing either 
of the diffusion differential equations 4, 5, or 9. 

3. Field Equations for Simulating Fluid Flow Over a Com- 
puter-Generated Wing Surface 

As discussed above, a fluid flow can be simulated using 
55 both inviscid and viscous CFD simulation modules to 
improve computing efficiency. The explanation below 
describes how the CFD simulation modules can be applied to 
different regions of fluid flow and provides exemplary field 
equations for simulating the dynamics of the fluid flow. 

60 FIGS, la and lb depict a cross-sectional representation of 
a fluid flow over a wing surface 702 classified by two regions: 
a free stream region 704 and a boundary-layer region 706. 

As shown in FIG. la, the free-stream region 704 may be 
located away from the wing surface 702. However, the free 
65 stream region 704 may be close to the wing surface 702 in 
areas where the boundary layer region 706 is thin or has yet to 
develop. See, for example, the portion of the fluid flow in FIG. 
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la near the leading edge of the wing surface 702. The free 
stream region 704 is usually characterized as having a rela- 
tively uniform velocity profile 712. When there is a uniform 
velocity profile 712, internal shear forces acting on a fluid 
may be relatively small, and therefore viscous contributions 
to the fluid dynamics can be minimized or ignored. 

As mentioned previously, a CFD simulation module that 
ignores viscous effects may also be called an inviscid CFD 
simulation module. Equation 13, below, provides an exem- 
plary field equation for an inviscid CFD simulation module. 
Equation 13, also called the Euler method, represents the 
conservation of mass, conservation of three components of 
momentum, and conservation of energy: 


6m 6f x 6f y 6f z Equation 13 

6r + 6x + 6y + 6z 

where 


and Design System for Low Reynolds Number Airfoils,” pp. 
1-12, Proceedings of the Conference on Low Reynolds Num- 
ber Aerodynamics (T. J. Mueller ed., Univ. of Notre Dame, 
Notre Dame, Ind., 1989). 

Equation 14, below, represents a boundary-layer integral 
momentum equation for compressible flow: 


to 


de 

dx 


, 6 du e 

+ (2 + h-m;)- — 

u e dx 


El 
2 ’ 


Equation 14 


where 0 is the momentum thickness, H is the shape factor, M e 
is the boundary-layer edge Mach number, u e is the boundary - 
layer edge velocity, and Cyis the skin friction coefficient. 

Equation 15, below, represents a boundary-layer kinetic 
energy integral equation: 
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+ (2H** + H*(l -H)) 


6 du e 
u e dx 


= 2 C D - H* - 


Equation 15 


25 As used in equations 14 and 1 5, above, shape factors H, H*, 
and FI** are defined as: 


30 


<5* e * 6 ** 

H = —H‘ = —H‘- = — ; 

6 9 9 


displacement thickness 6 * is defined as: 


where u, v, and w are components of the velocity vector, p is 
the pressure, p is the density, and E is the total energy per unit 35 
volume. Combining equation 13 with an equation of state 
(e.g., the ideal gas law), an inviscid CFD simulation module 
can predict the fluid properties for the free stream fluid region 
704. 

FIGS, la and lb also depict a boundary-layer region 706, 40 
located near a wing surface 702. The boundary layer region 
706 is typically characterized by a sharply increasing velocity 
profile 708. Skin friction causes the velocity of the fluid flow 
very close to the wing surface to be essentially zero, with 
respect to the surface. A sharply increasing velocity profile 45 
708 develops as the velocity increases from a near-zero veloc- 
ity to the free stream velocity. The sharply increasing velocity 
profile 708 in the boundary-layer region 706 creates shear 
forces within the boundary-layer region 706. Due to the inter- 
nal shear forces, viscous properties of the fluid influence fluid 50 
flow in the boundary-layer region 706. Therefore, a simula- 
tion of the fluid flow in the boundary-layer region 706 should 
account for viscous contributions to the flow dynamics. In 
some cases, the fluid flow in boundary-layer region 706 may 
be characterized as turbulent flow (region 710). Due to fluid 55 
vorticity, viscous properties of the fluid influence the fluid 
flow. Thus, a simulation of the turbulent flow should also 
account for viscous contributions to the flow dynamics. For 
purposes of this discussion, laminar and turbulent regions are 
treated as one boundary-layer region and simulated using a 60 
single CFD simulation module. 

As mentioned previously, a CFD simulation module that 
accounts for viscosity may also be called a viscous CFD 
simulation module or a boundary-layer CFD simulation mod- 
ule. Below, exemplary field equations for a boundary-layer 65 
CFD simulation module are provided according to a Drela 
boundary-layer technique. Drela, M. “XFOIL: An Analysis 


pv— y* 

Jo V Pe u e > 


momentum thickness 0 is defined as: 



kinetic energy thickness 0* is defined as: 


0 * = 


D‘~&; 




density thickness 5** is defined as: 


6- = 


r< 


1- - -Jy; 

pP U, 


skin friction coefficient C^is defined as: 
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and dissipation coefficient C D is defined as: 


Cn = 


if 


’ du 

T a~ dy - 

dy 


5 


Solving equations 14 and 15 for local velocity u and den- 
sity p, the boundary-layer CFD simulation module can pre- 
dict the fluid properties for portions of the fluid flow within 0 
the boundary-layer region 706. Additionally, characteristics 
of the boundary layer, including boundary-layer thickness, 
can also be determined once the fluid properties are known. 

In some cases, the boundary-layer CFD simulation module 
returns a transpiration flux value representative of a bound- 15 
ary-layer region having a predicted thickness 5*. A transpi- 
ration flux is a fictitious fluid flow into or out of the computer- 
generated aircraft surface. The transpiration flux displaces a 
portion of the free-flow region 704 in the same way as would 
the presence of a boundary layer having a predicted thickness 20 
6 *. This allows the boundary-layer region to be discretized 
using a set of boundary-layer prediction points instead of 
using small volumes or fluid cells. In general, as the magni- 
tude of the transpiration flux increases, the fictitious fluid flow 
increases, simulating a thicker boundary layer. In some cases, 2 5 
the transpiration flux can be used to create a fictitious flow of 
air into the aircraft surface (negative flux), thereby simulating 
a boundary layer having a reduced thickness. 

The transpiration flux can be determined using the output , f) 
of the Drela boundary-layer technique described in equations 
14 and 15, above. For example, the transpiration flow velocity 
W iw of the transpiration flux can be determined using: 


w iw 


1 d 
Piw <$S 


(piM w < n 


35 

Equation 16 


where p iw is the density of the fluid flow at the aircraft surface, 
U iw is the velocity of the fluid flow at the aircraft surface, 5* 40 
is the computed boundary-layer displacement thickness, and 
s is an arc length along the computer-generated aircraft sur- 
face. In some cases, the arc length s is the average distance 
between boundary-layer prediction points on the computer- 45 
generated aircraft surface. Equation 1 6 is taken from Lock, R. 

C., and Williams, B. R., “Viscous-Inviscid Interactions in 
ExtemalAerodynamics,”Prog.fferas/?ace&7.,Vol. 24, 1987, 
pp. 51-171. Thus, the transpiration mass flux (density p !VV 
times the transpiration flow velocity W !14 .) is equal to the rate 50 
of change of the product of the local density p !W , local velocity 
U, w , and boundary-layer displacement thickness 8* along the 
solution strip. A finite difference method can be used to com- 
pute the derivative in equation 16. For example, neighboring 
boundary-layer prediction points along the surface of the 55 
computer-generated aircraft surface can be used with a sec- 
ond order, backward Lagrange polynomial formulation to 
compute the derivative values. 

The results of the inviscid CFD simulation module (using, 
for example, equation 13) are used by a viscous or boundary 60 
layer CFD module (using, for example, equations 14, 15, and 
16) to predict the boundary-layer fluid properties for each 
boundary-layer prediction point. The boundary-layer fluid 
properties (e.g., boundary-layer thickness or transpiration 
flux) can then be used by the inviscid CFD simulation module 65 
to predict an updated or refined set of inviscid fluid property 
values. In this way, the inviscid and viscous CFD simulation 


modules can be used to iterate or refine the results of the 
computer-generated simulation. 

4. Exemplary Results 

FIGS. 8 through 10 depict exemplary pressure distribu- 
tions over a computer-generated wing surface. The following 
discussion compares results obtained using isotropic diffu- 
sion to results obtained using anisotropic diffusion. 

In general, isotropic diffusion may be adequate for an 
aircraft surface in flow conditions producing relatively homo- 
geneous boundary-layer fluid properties. As discussed above, 
this condition may be indicated by regions where the pressure 

gradient vector Vp is small or nearly zero. 

However, in situations where the fluid properties change 
rapidly over an aircraft surface (e.g., the surface of a wing), 
isotropic diffusion may blend or distort regions of the surface 
having different fluid properties. For example, a wing surface 
having a rapidly changing pressure distribution due to the 
formation of a shock wave, the shape of the pressure distri- 
bution may become distorted or blurred when using an iso- 
tropic diffusion technique. 

FIG. 8 depicts an exemplary pressure distribution pro- 
duced by diffusing a fluid pressure over a computer-generated 
wing surface using: an initial condition shown in FIG. 5; fixed 
boundary conditions for strip areas 404; and isotropic diffu- 
sion discretized with a finite-volume solver. The pressure 
distribution of FIG. 8 can be compared to an exemplary 
pressure distribution depicted in FIG. 9 representing the pres- 
sure distribution predicted using an inviscid or Euler-based 
CFD simulation module. 

As shown in FIG. 8, the pressure is lowest (designated by 
darker regions) in locations close to the strip areas where the 
fluid pressure was initially known. As explained above, a 
boundary-layer CFD simulation module can be used to deter- 
mine the boundary-layer fluid properties associated with the 
surface-mesh polygons in the strip areas. The uneven distri- 
bution along the length of the wing indicates that diffusion 
perpendicular to the strip areas 404 is not sufficient to accu- 
rately propagate the pressure fluid property. Increasing the 
relative diffusion rate along this direction would reduce this 
effect. Also, comparing FIG. 8 and FIG. 9, the pressure dis- 
tribution appears to be elongated along the chord of the wing 
in locations close to the strip areas 404. This indicates that 
there may be too much diffusion between different regions of 
pressure, resulting in distortion of the pressure regions. 
Decreasing the relative diffusion rate along this direction may 
reduce this effect. However, using the isotropic diffusion 
described in equation 9, above, the diffusion rate must be 
equal in all directions. 

Anisotropic diffusion provides a variable diffusion rate (or 
conduction coefficient) that depends on the direction of dif- 
fusion (see, for example, equations 4 and 5, above). Using an 
anisotropic diffusion technique, the diffusion rate can be tai- 
lored to allow more diffusion in directions where propagation 
is desired and less diffusion in other directions to minimize 
distortion of the fluid property distribution. In particular, it is 
desirable to limit diffusion in directions where the fluid prop- 
erties are changing most rapidly. In directions where the fluid 
properties are more homogenous, diffusion can be acceler- 
ated, allowing a more evenly-distributed propagation of the 
fluid property. 

FIG. 10 depicts an exemplary pressure distribution over a 
wing surface produced by propagating or interpolating the 
fluid pressure overthe wing surface using: an initial condition 
shown in FIG. 5; fixed boundary conditions on strip areas 
404; and pressure-gradient-driven anisotropic diffusion dis- 
cretized using a finite-volume solver. Note that the property 
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distribution is similar to the inviscid flow pressure distribu- 
tion shown in FIG. 9. In contrast to the property distribution 
shown in FIG. 8, the anisotropic diffusion does not distort the 
property distribution and produces a result that is more like 
the actual distribution over the exemplary aircraft surface. 5 
5. Computer and Computer Network System 

The embodiments described herein are typically imple- 
mented as computer software (computer-executable instruc- 
tions) executed on a processor of a computer system. FIG. II 
depicts an exemplary computer system 1100 configured to to 
perform any one of the above-described processes. Computer 
system 1100 may include the following hardware compo- 
nents: processor 1102, data input devices (e.g., keyboard, 
mouse, keypad) 1104 , data output devices (e.g., network con- 
nection, data cable) 1106 , and user display (e.g., display 15 
monitor) 1108 . The computer system also includes non-tran- 
sitory memory components including random access 
memory (RAM) 1110 , hard drive storage 1112 , and other 
computer-readable storage media 1114 . 

Processor 1102 is a computer processor capable of receiv- 20 
ing and executing computer-executable instructions for per- 
forming any of the processes described above. Computer 
system 1100 may include more than one processor for per- 
forming the processes. The computer-executable instructions 
may be stored on one or more types of non-transitory storage 25 
media including RAM 1110 , hard drive storage 1112 , or other 
computer-readable storage media 1114 . Other computer- 
readable storage media 1114 include, for example, CD-ROM, 
DVD, magnetic tape storage, magnetic disk storage, solid- 
state storage, and the like. 30 

FIG. 12 depicts an exemplary computer network for dis- 
tributing the processes described above to multiple computers 
at remote locations. One or more servers 1210 may be used to 
perform portions of the process described above. For 
example, one or more servers 1210 may store and execute 35 
computer-executable instructions for receiving information 
for generating a computer-generated simulation. The one or 
more servers 1210 are specially adapted computer systems 
that are able to receive input from multiple users in accor- 
dance with a web-based interface. The one or more servers 40 
1210 are able to communicate directly with one another using 
a computer network 1220 including a local area network 
(LAN) or a wide area network (WAN), such as the Internet. 

One or more client computer systems 1240 provide an 
interface to one or more system users. The client computer 45 
systems 1240 are capable of communicating with the one or 
more servers 1210 over the computer network 1220 . In some 
embodiments, the client computer systems 1240 are capable 
of running a web browser that interfaces with a Web-enabled 
system running on one or more server machines 1210 . The 50 
Web browser is used for accepting input data from the user 
and presenting a display to the user in accordance with the 
exemplary diffusion techniques described above. The client 
computer 1240 includes a computer monitor or other display 
device for presenting information to the user. Typically, the 55 
client computer 1240 is a computer system in accordance 
with the computer system 1100 depicted in FIG. 11 . 

Although the invention has been described in considerable 
detail with reference to certain embodiments thereof, other 
embodiments are possible, as will be understood by those 60 
skilled in the art. 

We claim: 

1. A computer-implemented method of generating a fluid- 
flow simulation over a computer-generated aircraft surface, 65 
the surface comprised of a surface mesh of polygons, the 
method comprising: 
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obtaining a boundary-layer fluid property for a subset of 
the polygons of the surface mesh of polygons, the subset 
of polygons representing a portion of the computer- 
generated aircraft surface; 

determining, using a computer, a pressure-gradient vector 
for a selected polygon, the selected polygon belonging 
to the surface mesh of polygons but not one of the subset 
of polygons for which a boundary-layer fluid property 
has been obtained; 

defining a minimum diffusion rate along the direction of 
the pressure-gradient vecto of the selected polygon; 

defining a maximum diffusion rate in a direction perpen- 
dicular to the direction of the pressure-gradient vector of 
the selected polygon; 

defining a diffusion-path vector from a point in the selected 
polygon to a neighboring point in a neighboring poly- 
gon, the neighboring polygon belonging to the subset of 
polygons for which a boundary -layer fluid property has 
been obtained; and 

determining, using the computer, an updated fluid property 
for the selected polygon using a variable diffusion rate, 
the variable diffusion rate based on the minimum diffu- 
sion rate, maximum diffusion rate, and angular differ- 
ence between the diffusion-path vector and the pressure- 
gradient vector. 

2. The computer-implemented method of claim 1, wherein 
the variable diffusion rate is equal to the square root of the 
sum of the square of the minimum diffusion rate times cosine 
theta and the square of the maximum diffusion rate times sine 
theta, wherein theta is the angle between the direction of the 
diffusion-path vector and the direction of the pressure-gradi- 
ent vector. 

3. The computer-implemented method of claim 1, wherein 
determining the updated fluid property for the selected poly- 
gon is performed using an explicit discretization comprising: 

determining a distance between the point in the selected 
polygon and the neighboring point in the neighboring 
polygon; 

determining a property difference between an initial fluid 
property of the selected polygon and a boundary-layer 
fluid property of the neighboring polygon; and 

determining an updated fluid property by adding the initial 
fluid property value to the variable diffusion rate multi- 
plied by the property difference divided by the distance. 

4 . The computer-implemented method of claim 1 , wherein 
determining the updated fluid property for the selected poly- 
gon is performed using an implicit discretization comprising: 

obtaining coefficient values for a coefficient matrix, the 
coefficient values based on discretized spatial deriva- 
tives associated with the surface mesh of polygons; 

determining a set of linear equations using the coefficient 
matrix, wherein each linear equation of the set of linear 
equations corresponds to a surface-mesh polygon of the 
surface mesh of polygons; 

determining an updated property value associated with the 
selected polygon by simultaneously solving the set of 
linear equations, the updated property value represent- 
ing the boundary-layer fluid property of the selected 
polygon. 

5. The computer-implemented method of claim 1, wherein 
the point in the selected polygon is a centroid of the selected 
polygon and the neighboring point in the neighboring poly- 
gon is a centroid of the neighboring polygon. 

6. The computer-implemented method of claim 1, wherein 
the boundary-layer fluid property is a boundary-layer dis- 
placement thickness. 
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7. The computer-implemented method of claim 1, wherein 
the boundary-layer fluid property is a transpiration flow vec- 
tor, the transpiration flow vector having a direction and a rate 
of fluid flow originating from the aircraft surface so as to 
represent the presence of a boundary-layer fluid flow. 

8. The computer-implemented method of claim 1, wherein 
obtaining the boundary-layer fluid property comprises: 

determining an inviscid fluid property for a set of fluid cells 
corresponding to the subset of the polygons using an 
inviscid fluid simulation that does not simulate fluid 
viscous effects; determining a boundary-layer predic- 
tion point for each polygon of the subset of the polygons; 
and 

determining a boundary-layer fluid property for each 
boundary-layer prediction point using inviscid flow data 
and a boundary-layer simulation that simulates fluid 
viscous effects. 

9. The computer-implemented method of claim 1, wherein 
the aircraft surface is the surface of a wing. 

10 . A computer-implemented method of generating a fluid- 
flow simulation over a computer-generated aircraft surface, 
the surface comprised of a surface mesh of polygons, the 
method comprising: 

determining, using a computer, a pressure-gradient vector 
for a selected polygon, the selected polygon belonging 
to the surface mesh of polygons; 

defining a minimum diffusion rate along the direction of 
the pressure-gradient vector of the selected polygon; 

defining a maximum diffusion rate in a direction perpen- 
dicular to the direction of the pressure-gradient vector of 
the selected polygon; 

defining a diffusion-path vector from a point in the selected 
polygon to a neighboring point in a neighboring poly- 
gon, the neighboring polygon belonging to the surface 
mesh of polygons; and 

determining, using the computer, an updated fluid property 
for the selected polygon using a variable diffusion rate, 
the variable diffusion rate based on the minimum diffu- 
sion rate, maximum diffusion rate, and angular differ- 
ence between the diffusion-path vector and the pressure- 
gradient vector. 

11. A non-transitory computer-readable medium storing 
computer-readable instructions for generating a fluid-flow 
simulation over a computer-generated aircraft surface, the 
surface comprised of a surface mesh of polygons, the instruc- 
tions performing a process comprising: 

obtaining a boundary-layer fluid property for a subset of 
the polygons of the surface mesh of polygons, the subset 
of polygons representing a portion of the computer- 
generated aircraft surface; 

determining a pressure-gradient vector for a selected poly- 
gon, the selected polygon belonging to the surface mesh 
of polygons but not one of the subset of polygons for 
which a boundary-layer fluid property has been 
obtained; 

defining a mini mum diffusion rate along the direction of 
the pressure-gradient vector of the selected polygon; 

defining a maximum diffusion rate in a direction perpen- 
dicular to the direction of the pressure-gradient vector of 
the selected polygon; 

defining a diffusion-path vector from a point in the selected 
polygon to a neighboring point in a neighboring poly- 
gon, the neighboring polygon belonging to the subset of 
polygons for which a boundary-layer fluid property has 
been obtained; and 

determining an updated fluid property for the selected 
polygon using a variable diffusion rate, the variable dif- 
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fusion rate based on the minimum diffusion rate, maxi- 
mum diffusion rate, and angular difference between the 
diffusion-path vector and the pressure-gradient vector. 

12. The computer-readable medium of claim 11, wherein 
5 the variable diffusion rate is equal to the square root of the 

sum of the square of the minimum diffusion rate times cosine 
theta and the square of the maximum diffusion rate times sine 
theta, wherein theta is the angle between the direction of the 
diffusion-path vector and the direction of the pressure-gradi- 
10 ent vector. 

13. The computer-readable medium of claim 11, wherein 
determining the updated fluid property for the selected poly- 
gon is performed using an explicit discretization comprising: 

determining a distance between the point in the selected 
15 polygon and the neighboring point in the neighboring 
polygon; 

determining a property difference between an initial fluid 
property of the selected polygon and a boundary-layer 
fluid property of the neighboring polygon; and 
20 determining an updated fluid property by adding the initial 
fluid property value to the variable diffusion rate multi- 
plied by the property difference divided by the distance. 

14. The computer-readable medium of claim 11, wherein 
determining the updated fluid property for the selected poly- 

25 gon is perforated using an implicit discretization comprising: 

obtaining coefficient values for a coefficient matrix, the 
coefficient values based on discretized spatial deriva- 
tives associated with the surface mesh of polygons; 

determining a set of linear equations using the coefficient 
30 matrix, wherein each linear equation of the set of linear 
equations corresponds to a surface-mesh polygon of the 
surface mesh of polygons; 

determining an updated property value associated with the 
selected polygon by simultaneously solving the set of 
35 linear equations, the updated property value represent- 
ing the boundary-layer fluid property of the selected 
polygon. 

15. The computer-readable medium of claim 11, wherein 
the point in the selected polygon is a centroid of the selected 

40 polygon and the neighboring point in the neighboring poly- 
gon is a centroid of the neighboring polygon. 

16. The computer-readable medium of claim 11, wherein 
the boundary-layer fluid property is a boundary-layer thick- 
ness value, the boundary-layer thickness value representing a 

45 distance from the aircraft surface where fluid viscous effects 
can be ignored. 

17. The computer-readable medium of claim 11, wherein 
the boundary -layer fluid property is a transpiration flow vec- 
tor, the transpiration flow vector having a direction and a rate 

50 of fluid flow originating from the aircraft surface so as to 
represent the presence of a boundary-layer fluid flow. 

18. The computer-readable medium of claim 11, wherein 
obtaining the boundary-layer fluid property comprises: 

determining an inviscid fluid property for a set of fluid cells 
55 corresponding to the subset of the polygons using an 
inviscid fluid simulation that does not simulate fluid 
viscous effects; determining a boundary-layer predic- 
tion point for each polygon of the subset of the polygons ; 
and 

60 determining a boundary-layer fluid property for each 
boundary-layer prediction point using inviscid flow data 
and a boundary-layer simulation that simulates fluid 
viscous effects. 

19. The computer-readable medium of claim 11, wherein 
65 the aircraft surface is the surface of a wing. 

20. A non-transitory computer-readable medium storing 
computer-readable instructions for generating a fluid-flow 
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simulation over a computer-generated aircraft surface, the 
surface comprised of a surface mesh of polygons, the instruc- 
tions performing a process comprising: 

determining a pressure-gradient vector for a selected poly- 
gon, the selected polygon belonging to the surface mesh 5 
of polygons; 

defining a minimum diffusion rate along the direction of 
the pressure-gradient vector of the selected polygon; 
defining a maximum diffusion rate in a direction perpen- 
dicular to the direction of the pressure-gradient vector of to 
the selected polygon; 

defining a diffusion-path vector from a point in the selected 
polygon to a neighboring point in a neighboring poly- 
gon, the neighboring polygon belonging to the surface 
mesh of polygons; and 15 

determining an updated fluid property for the selected 
polygon using a variable diffusion rate, the variable dif- 
fusion rate based on the minimum diffusion rate, maxi- 
mum diffusion rate, and angular difference between the 
diffusion-path vector and the pressure-gradient vector. 20 



